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We analyze the fluctuations of an electronic heat current 
across an idealized molecular junction. Apart from an 
induced, finite inter-electrode heat current the focus will 
be on the spectral features of the resulting heat current 
fluctuations. By use of the Green function method we de- 
rive an explicit expression for the power spectral density 
of the emerging heat current fluctuations. The complex 
expression simplifies considerably in the limit of zero 
frequency, yielding the noise intensity of the heat cur- 
rent. This spectral density for the heat fluctuations still 
depends on the frequency in the zero-temperature limit, 
assuming different asymptotic behaviors in the low- and 
high-frequency regions. We further address some sub- 
tleties and open problems from an experimental view- 
point. 



1 Introduction The experimental activities over the 
last fifteen years in investigating transport across molecu- 
lar junctions UJ2) has triggered a wave of intense research 
in both communities, theoretically QUEM) an d as we U 
as in experiment |7 <X 9 l()|. Single molecule electronics is 
still considered as a possible candidate for the substitution 
of the silicon-based elements in the information process- 
ing technology 08). This fact makes studies of molecular 
transport properties very appealing. Likewise, such molec- 
ular junctions seemingly have advantages in the context 
of energy-related applications. This is due to their poten- 
tial for hybrid solid-state molecular structure which yields 
novel interface properties and the abundant selection of 




Gate 




(color online) Sketch of a molecular junction which transfers 
thermal energy carried by electrons from a hot electrode TL 
across an idealized molecular junction towards a neighboring 
cold one Tr. Here, the inter-electrode electronic level eo can 
be controlled continuously. 
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possible molecules and electrode materials which allow to 
tailor specific properties on demand. In particular the issue 
of thermoelectric conversion and photovoltaic conversion 
processes continues to attract interest in pursuing research 
in molecular electronics. 

Apart from the current-voltage characteristics [Q]|9] 
[TO), it is also possible, for example by use of full count- 
ing statistics 111 llll2|[T3l[T4l . to extract information about 
the fluctuations of the electric current flowing through a 
molecular wire II15II16|[T7II18II191 . In this context, the issue 
of heat or phonon transport across such molecular junc- 
tions plays a non-negligible role which may well crucially 
impact the electronic transport features, including the over- 



Copyright line will be provided by the publisher 



2 



First author et al.: Short title 



all stability issue of the molecular junction setup. It must 
also be realized that the experimental exploration of heat 
transport presents not at all a straightforward task. With 
this work we shall mainly explore the fluctuations of heat 
current caused by the transferring electrons. 

With the systems of interest operating on the nanoscale 
the corresponding heat current fluctuations can become 
sizable. This may be so even in situations where the av- 
erage heat current is vanishing identically, as it is the 
case in thermal equilibrium with both interconnecting 
electrodes held at the same temperature. Moreover, the 
properties of noise correlation features, or likewise, its 
frequency-dependent spectral properties and, as well, its 
zero-frequency power spectrum, are in no obvious manner 
related to the mean value of the heat flow itself. Our goal is 
to obtain analytical estimates for its power spectral density 
(PSD), even at the expense that these may mainly apply to 
idealized setups only. With such a restriction these analyt- 
ical results may nevertheless be useful to appraise the role 
of heat current noise in more realistic molecular junctions. 
It is further of interest to have an estimate available when 
devising molecular circuitry for more complex tasks. 

Heat transport across a molecular structure which links 
two electrodes is induced by a difference of the corre- 
sponding lead temperatures, i.e. Tl at the left-sided elec- 
trode and Tr at the right-sided electrode. The physics of 
heat transfer generally involves both electrons and phonons 
and their mutual interaction [6. 20. 21. 22. 23. 24. 25. 26 ,27 , 
28,29]. Therefore, the intensity of heat flow carried across 
the wire should be addressed with care, with the need to 
distinguish between heat transfer mediated either by elec- 
trons or phonons, or a combination of both. If phonons are 
mainly at work this topic relates to the new field of phonon- 
ics l30l . a novel research area which may lead to new cir- 
cuit elements, such as molecular thermal diodes, thermal 
transistors, thermal logic gates, to name but a few H301I3U 
[32ll33U34U35lF36l . Here also, the size of fluctuations in heat 
current does matter; this is so because those may well turn 
out to be deleterious to intended information processing 
tasks. 

Energy transport mediated by electrons which relates 
at the same time to charge transfer: electrons are moving 
from lead to lead, carrying not only charge but also en- 
ergy. However, the amount of energy transferred by a sin- 
gle electron, unlike to its charge, is not quantized II37II381 . 
In contrast to the studies that examine the average heat flow 
much less attention, however, has been paid to the issue and 
impact of fluctuations of the accompanying heat flow. In a 
prior work [39 ] the heat transport through a ballistic quan- 
tum wire has been considered in the Luttinger-liquid limit, 
by neglecting the discreteness of the wire's energy spec- 
trum. Likewise, with Refs. H40II411 . the PSD of the heat 
current fluctuations has been derived within the scattering 
theory, under the assumption that the electrons are trans- 
mitted (reflected) at the same rates, independently of their 
actual energies. The results of the last two papers, how- 



ever, are distinctive because they have shown that the noise 
characteristics of heat current exhibits a well-pronounced 
frequency dependence even in the zero-temperature limit. 

2 Molecular junction setup 

Here we shall consider the electronic heat current that pro- 
ceeds across a molecular wire composed of a single en- 
ergy level with the two leads held at constant tempera- 
tures. A preliminary short report of electronic heat noise 
features has been presented by us with Ref. J42]- Here we 
extend this study and present further details on the theoret- 
ical derivation of the noise expression and discuss the fea- 
tures of the noise of the emerging heat current over much 
broader parameter regimes and frequency regimes away 
from the zero-frequency result. In order to obtain analytic al 
tractable expressions we shall neglect electron-phonon in- 
teractions and, as well, electron-electron interactions. Such 
a simplification can be justified for tailored situations that 
involve a very short wire only. Then, the Coulomb interac- 
tion via a double occupancy shifts the energy far above the 
Fermi level so that its role in thermal transport can be ne- 
glected. Likewise, the electron dwell time is short as com- 
pared to the electron-phonon relaxation time scale. Note 
however, that in contrast to previous works [40j|4T] we ac- 
count here for the dependence of the transmission coeffi- 
cient on its electron energies, and, within the Green func- 
tion approach [5 29], derive an explicit expression for the 
PSD of the heat current fluctuations, S h (u>). In particular 
we show below that the net noise features of the heat cur- 
rent are quite distinct from their electronic counterpart. 

Our model molecular junction setup is depicted with 
Fig-E It is described by a Hamiltonian 

H = -Hwire + Pleads + ^contacts ■ (1) 

It contains three different contributions, namely the wire 
Hamiltonian, the leads and the wire-lead coupling, respec- 
tively. We consider here the regime of coherent quantum 
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Figure 1 (color online) Idealized setup of a molecular 
junction used in text: Two metal leads, filled with elec- 
tron gas, are connected by a single orbital Eq. The coupling 
strengths are determined by the constants Tl/r. The left 
lead is prepared at a higher temperature as compared to the 
opposite right lead, i.e. Tl > Tr. The chemical potential, 
/i, is the same for both leads so that no electric current due 
to a voltage bias is present. 
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transport whereby neglecting dissipation inside the wire. 
The wire is composed of a single orbital; i.e., 



sod'd , 



(2) 



at an energy Eq, with the fermionic creation and annihi- 
lation operators, and d. The energy level £o can be 
tuned by applying a gate voltage. This so set up idealized 
model allows for explicit analytical calculations. Physi- 
cally, it mimics a double barrier resonant tunneling struc- 
ture GaAs/AlajGai-^-structure of the type considered for 
electronic shot noise calculations in Ref. B3l . herein trun- 
cated to a single Landau level. As commonly done, the 
electrodes are modeled by reservoirs, composed of ideal 
electron gases, i.e., 



-Pleads — 2_j Se< ! C eq C£ 1 ' 



(3) 



where the operator c\ q (ci q ) creates (annihilates) an elec- 
tron with momentum q in the £ ~L (left) or I =R (right) 
lead. We assume that the electron distributions in the leads 
are described by the grand canonical ensembles at the tem- 
peratures Tl/r ar, d with chemical potentials /Jl/r- With 
such ideal electron reservoirs we then obtain 



where 



(c\ c eg >) = 5u'5 qq 'ft{ei q ) , 



f t (e £q ) = U^-^)AbT, + 1 



(4) 



(5) 



denotes the Fermi function. 

We impose a finite temperature difference AT = XL — 
Tr and use identical chemical potentials, /il = /Ur = M f° r 
the electrodes. When an electron tunnels out from a lead, 
the energy E is transferred into the wire which presents 
the heat transfer, 8Q. Observing the value for the chemical 
potential, /i, it reads SQ = E — /.t. In the following we 
use that all the electron energies are measured from the 
chemical potential value /i, being set at fx = 0. 

The Hamiltonian which describes the tunneling events 
reads: 

-^contacts = ^ Vi q cj, q d + tl.C. . (6) 

This part mediates the coupling between the wire and the 
electrodes. Here, the notation h.c. denotes Hermitian con- 
jugate. The quantity Vi q is the tunnelling matrix element, 
and the tunneling coupling is characterized in general by a 
spectral density, 



r e {E) = 2TrY,\Vt q \ 2 S(E-e eq ) 



(7) 



In the following, we shall use a wide-band limit of the 
electrode conduction bands, setting Tg(E) := Tg. 



3 Power spectral density of electronic heat cur- 
rent fluctuations 

Working within the Heisenberg description of operators we 
present the derivation of the electronic heat current and the 
PSD of its corresponding fluctuations. The heat current is 
given by 



At 



(8) 



With our choice of chemical potentials, we find that 
the heat transfer operator is SQ(t) = E^, with the energy 
operator given by 



E, 



Lg c Lg c L<3 • 



(9) 



Its time derivative thus yields the operator for the heat cur- 
rent, reading: 



Ji{t) = -Y. 2j w L MVuci q w{t)\ 



(10) 



The heat current is positive valued when heat transport pro- 
ceeds from the hot left lead, i.e. Tl > Tr to the adjacent 
cold lead, see in Fig. 1. Deriving the above expression, we 
have employed the Heisenberg representation for the lead 
electron operators. The average current is obtained by the 
ensemble average (J^(t)). Because we consider only the 
heat transport by electrons, we have (J^(t)) = ^(J^(t)), 
due to the conservation of energy. We henceforth focus on 
the quantities derived with regard to the left lead. 

The quantum correlation function of heat current fluc- 
tuations is described by the symmetrized autocorrelation 
function, i.e., 

S h (t,t') = ±([AJ*(t),Aj£ (f')]+> , (11) 

with respect to the operator of the heat current fluctuation 

AJt(t) - j£(t) - (j£(t)) . (12) 

The heat current noise is described here with r = t — t' 
by the symmetrized quantum autocorrelation function 



S h (r)^l/2([Ajf(r),A4(0)] + ), 



(13) 



of the heat current fluctuation operator Ajf(s) = Jf(s) — 
(Jf(s)), where the anti-commutator [A, B] + = AB + BA 
ensures the hermitian property. 

With this work we throughout consider the asymptotic 
long time limit t — > oo when all transients are decayed. In 
this asymptotic limit the average heat current is stationary 
and the autocorrelation function of the heat current fluctua- 
tions becomes time-homogeneous; i.e. it is independent of 
initial preparation effects. It thus depends on the time dif- 
ference t = t — t' only. The Fourier transform yields the 
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power spectral density (PSD) S h (uj) for the heat current 
noise, i.e., 



S h (io) = S h (-io) = / dTe l " T S h (T) > 



(14) 



S h (Lu) is an even function in frequency and strictly semi- 
positive, in accordance with the Wiener- Khintchin theorem 
B4l . In the following we address positive values of the fre- 
quency, to > 0, only. 

The annihilation operators of the electrode states sat- 
isfy the Heisenberg equations of motion; i.e., 

C( q (t) = --e lq c lq (t) - -Vi q d{t) , (15) 
yielding the solution 

f dt'e-^-^dit') . (16) 

" J t 

Here, the first term on the right hand side describes the dy- 
namics of the free electrons in the leads, while the second 
term accounts for the influence of the molecule. 

The Heisenberg equation of the molecular annihilation 
operator is given by 

d(t) = - l -e Q d{t) - - £ V e * q c eq (t). (17) 

£q 

Upon inserting Eq. ( [TBI l into Eq. (fTTb . we obtain 

r L + r R 



d = -s d(t) - 



2h 



-d(i) + &(*) + &(*), (18) 



where we have defined the noise operator 



— ^££q(t — to) 



c iq (t ). (19) 



In addition, we have employed the definition ^} and used 
the wide-band limit. 

The noise quantity defined in Eq. ( fT9] l denotes operator- 
valued Gaussian noise, which is characterized by its mean 
and correlation properties, reading 



&(*)> = o 



(20) 



<&(0&(*)> = fc 



de 
2-kU' 



(21) 



This noise accounts for the influence of the states stem- 
ming from the electrodes I = L, R. 

Now the central problem is to solve the inhomogeneous 
differential equation (11 8t . Once we obtain the solution of 



Eq. (fT8l . we obtain also the solution for Eq. dT6b . the heat 
current (fTOt and also the power spectral density in Eq. ( \\A\ . 

To obtain the solution we follow the Green function 
approach in Ref. J5) and start with solving the following 
differential equation 

%£o ' rh + FR )G(t-t>) = 5(t-t>), (22) 



^dt + h 



2h 



followed by the application of the convolution d(t) = 
fG(t- t')(£ h (t') + £, R (t'))dt'. The solution of Eq. (fS| is 
thus given by: 

G{t) = 0{t) e - lEot / h -^ +r ^ t / 2h . (23) 
Then, the molecular operator in Eq. ([T8~t assumes the form 

^ q e eq ~ e + i{r L + r K ) 2 



In what follows we address solely the asymptotic prop- 
erties which are reached with the initial time of prepara- 
tion to — > — oo. This implies that average currents assume 
stationary values and correlation functions become time- 
homogeneous. With this expression and its Hermitian con- 
jugate, we obtain the occupation value of the molecular en- 
ergy level Eq as 

n E0 ={d^t)d(t)} 

V £q eyLp[iei q (t - t )/K\ 



E 



,7^ [e lq -e -i{r^+r K )/2] 
VI, , exp[-i£ e > q >(t - t )/h] 



{e£> q > — Eq 



+ i(r L + r R )/2] ( ^ ( * o)c ^' (io)) 



E 



\Viq\ 2 fe(seq) 



r ) 2 + (r L + r R )V4 ' 



(25) 



where we have employed the ensemble average, Eq. (0). 
We find that this occupation is determined first by the 
Fermi function of the leads weighted by the tunneling ma- 
trix elements Ve q , the difference among the energy levels 
of the lead states and between lead states and the molecular 
energy level £o. This occupation value is time-independent 
since there is no time-dependent external fields present. 

Upon substituting the result in Eq. (|24| | into Eq. H6i . 
we find for the operators in the electrodes 

CtS) = ci q (t )e-^ (t - to)/h 



E 



ce>q>{to) 



evg, -e +i(r h + r R )/2 



B[e Vl 



£ £q\ 



(26) 
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Figure 2 (color online) Average electronic currents (top 
row) for charge J el and heat flow J h and the zero- 
frequency values of corresponding PSDs of the accompa- 
nying current fluctuations (bottom row) for charge (left col - 
umn) and heat (right column) transport through the single- 
orbital wire as functions of orbital energy for Tr = Tl = 
r = 0.1 meV. The remaining parameters are Tl = 5.2 K, 
Tr = 3.2 K (solid lines). For equal chemical potentials 
/il = /iR = the heat flow in panel (b) vanishes for equal 
temperatures Tl = Tr; its fluctuations in equilibrium at 
the temperature Tl = Tr = 4.2 K are depicted versus 
orbital energy eo by the dashed line in panel (d). Figure 
adapted from Ref. J42). 



where, 



B(E) = V 



- iir6(E) , 



(27) 



and V denotes the integral principal value. In going from 
Eq. ( 124b to Eq. d26b we have used Sokhotsky's formula 
which states that rim e ->.o l/(x + ie) = V{l/x) — iiz5{x), 

where V(l/x) = f_ dx/x + dx/x, see in Ref. B31I . 



Next we insert Eq. d24l i and Eq. (|26| | into the heat cur- 
rent operator, Eq. (Hot , and by consequently taking the en- 
semble average, we obtain a Landauer-like formula for the 
heat current 028]|48] ; reading, 

(jh(t)> := Jh = i / dEEr (E)[h(E) - ME)}, 

(28) 

where the transmission coefficient 



T(E) = r h r R /[(E - e Q ) 2 + r 2 } 



(29) 



is energy-dependent. 

The expression for the electric Seebeck current [0 
reads very similar to Eq. (f28), except for its absence of the 
energy multiplier E in the integral in the rhs of Eq. (l28l . 
This seemingly small difference changes, however, the 
physics of the transport through the wire, because the mul- 
tiplier inverts the symmetry of the integral. Namely, the 
Seebeck current is an antisymmetric function of orbital en- 
ergy and vanishes when the orbital energy level is aligned 
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Figure 3 Power spectral density of the heat current noise 
as a function of the frequency ui at temperatures Tl = 
6K,Tr = 2K. The other parameters are eo = and 
r = 0.1 meV. 



to the chemical potentials of the leads, Fig. |2{a), while 
the heat current is a symmetric function and acquires a 
nonzero value at eo = 0, Fig. |2f b). 

3.1 Explicit results and discussion 

Upon combining Eq. (fl4] i and Eq. ([Tol l, we end up after a 
cumbersome evaluation with the nontrivial expression for 
the PSD of electronic heat current noise. Due to the com- 
plexity of this resulting expression the physics it inherits 
is not very illuminative. Nevertheless, we depict it here as 
given in our preliminary report lH2l . reading: 




E±-) T(E)T(E±f2) 



r 2 [E(E - e ) -(E± n)(E ± fi — e )Y 
[{E — eo) 2 + r 2 } [(E ±o — e ) 2 + r 2 } 

x h(E)J L (E ± Q) 

E±jj T(E)T(E ± f2)f R (E)J B (E ± 

r 2 T{E ± Q) 
(E-s ) 2 + r 2 

I- E 2 K(E)T(E ±f2)T ^ES1T(E)T(E ± SI) 
x h(E)J R (E ± Q) 

a 









(4) 







(E±f2) [±-j T{E)T{E±Q) 
(E ± ST) 2 K(E ± Q)T(E) 

n\ r 2 T(E±n) 



E±° 
2 



T 2 



(E - e ) 2 + T 2 



x f R (E)f L (E±f2) , 



(30) 
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Figure 4 Power spectral density of heat current noise at frequency uj = 2.16 x 10 13 Hz (which is the Debye cut-off 
frequency of gold) as a function of temperature difference with (a) weak molecule-wire coupling T = 0.1 meV or (b) 
strong molecule-wire coupling T = 10 meV. The other employed parameters are Tr = 300 K and Eq = 0. 



wherein we abbreviated Q = fiui, f = 1 — / with 7Z(E) 
I — T(E) denoting the reflection coefficient. Below we 
consider the case of symmetric coupling between the wire 
and the leads, Tl = Tr = T. It should also be men- 
tioned that its explicit verification is experimentally not 
straightforward. Typically such a PSD can be measured 
only indirectly via a single-time measurement of a tai- 
lored linear response function via a corresponding, gen- 
erally non-equilibrium quantum fluctuation-dissipation re- 
lation, which connects this response function with a cor- 
responding quantum two-time correlation expression lH4l . 
Put differently, this tailored response function is required to 
relate precisely to our quantum correlation of heat current 
fluctuations in d 1 3t . This is so because an explicit quantum 
measurement at two different times typically will impact 
(i.e. it will generally alter) the purely theoretically deter- 
mined quantum correlation expression in Eq. (TBI , for fur- 
ther details see in Refs. H46II47I . 

In Figure [3] we depict the dependence of the PSD of 
heat current noise versus frequency at the temperatures 
T L = 6 K, T R = 2 K. We find that the PSD depicts differ- 
ent power laws in different frequency regions and increases 
monotonically. 

Moreover, we find that the spectral density strength 
T of the wire-lead coupling can change the dependence 
of the PSD of heat current noise on the parameters. In 
Figure |4] we depict the PSD of heat noise as a function of 
the temperature difference AT in the case of weak and 
strong wire-lead couplings. With weak coupling, the PSD 
is smaller by one order of magnitude and only slightly 
depends on AT. In contrast, it increases very fast with AT 
when the coupling is very strong. According to Eq. ( 129) , 
the transmission coefficient becomes wider when T is 
larger, such that more electrons, whose energies deviate 
stronger from the chemical potential, are allowed to trans- 
port across the molecular junction. Therefore, the PSD 
becomes strongly enhanced and depends sensitively on 



AT. 

3.2 Zero frequency noise power 

The PSD of heat current noise at zero frequency ui = 
simplifies considerably, assuming an appealing form 

5 h (^ = 0;T L ,T R ) 

~ J dEE 2 [T(E)[ME)[l-fr(E)} 

+ f R (E)[l-f R (E)}) 

+ T(T)[1-T(T)][/ L (T)-/ R (T)] 2 ] . (31) 

Let us also contrast this result with the zero-frequency PSD 
of the fluctuations displayed by the nonlinear, accompany- 
ing Seebeck electric current. The latter reads EE]: 

5 cl ( W = 0;T L ,T R ) 
~ J dE[T(E)[f L (E)[l-f L (E)} 
+ f R (E)[l - f R (E)}) 

+ T(E)[1 -T(E)][f L (E) - f R (E)] 2 } , (32) 

Most importantly, the zero-frequency PSD for heat cur- 
rent in Eq. (l3~TT l differs by the energy factor E 2 within the 
integrand. Although this distinction seemingly appears mi- 
nor and may even be guessed beforehand without going 
through the laborious task of doing a theoretical rigorous 
derivation from which this limit derives from the final re- 
sult in Eq. (130b . It must be emphasized nevertheless that the 
two experessions lead to tangible differences. Particularly 
note the different behavior of the electronic and heat noise 
PSDs versus the tunable energy level £o as depicted with 
Fig. [2] and in Fig. [5] While the zero-frequency component 
of the electric PSD at oj = exhibits a maximum at Eq = 
its heat current PSD possesses instead a local minimum at 
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Figure 5 (color online) Power spectral density of the heat current noise at zero frequency u> = 0, (left panel) and power 
spectral density of the electric current noise (right panel) as functions of the wire orbital site energy Eq and wire-lead 
coupling strength T. The parameters employed are Tl = G-2K and Tr = "2.2K. 



this value. These two PSDs for charge current and heat cur- 
rent are compared in Fig. |5]over wide regimes of electronic 
orbital energy £o an d lead-molecule strength T. 

These differences originate from the salient feature 
that the two transport mechanisms for charge and the en- 
ergy are different. The electric current is quantized by the 
electron charge, e, while, in contrast, the energy carried 
by the electron is continuous and can assume principally 
an arbitrary value. Notably, the main contribution to the 
electronic noise power across the wire stems from those 
electrons occupying energy levels around the chemical po- 
tential fj, = 0. When £o deviates from the chemical poten- 
tial, increasingly less electrons participate in the transport. 
The flow of electron becomes diminished, and since both, 
the electric current and the electric noise are insensitive 
to the electron kinetic energies, they both decrease with 
increasing |erj|. This scenario differs for heat flow: There, 
the deviation from the chemical potential increases the 
possibility that successive electrons will carry different en- 
ergies. This in turn causes an increase of heat current noise. 
With further deviation of the orbital energy from the chem- 
ical potential, the occupancy difference [/l(-E') — fn(E)] 
decreases monotonically; consequently the heat flow noise 
power decreases again. 

3.3 Electronic heat current noise in thermal equi- 
librium 

Finally, let us focus on the thermal equilibrium heat cur- 
rent noise properties which are attained when the two tem- 
peratures are set equal, i.e. if Xl = Tr. In this case the 
average heat current vanishes identically, while its fluctu- 
ations remain finite. The zero-frequency spectra of both 
noise spectra for heat and electric current noise increase 
upon increasing the coupling strength T. This is so be- 
cause the transmission probability increases. The corre- 



sponding noise intensities assume nonzero values in equi- 
librium, however, as depicted with Fig. [2] see in panel 2(d). 

The properties at vanishing temperature, Tl = Tr, = 
0, are even more subtle. Here, the heat current PSD still 
depends on frequency. This dependence originates from 
quantum fluctuations where virtual transitions of electrons 
directly from lead-to-lead occur H401I411 . The Fermi dis- 
tribution equals the Heaviside step function in this case. 
Therefore, the contributions to the integrand in Eq. [30] 
stems from the interval \—D, 0]. After an integration of 
Eq. d30l l. one finds for the frequency dependent PSD the 
expression: 



0) 



S h (u,T L = T R 

r ( ( n 

[(2Qf - 2r 2 ] arctan 



+ 2 or 



r 



1 +loe 



r 4 



(n 2 + n) 



fi = hui. 



(33) 



In the limit r — > oo the zero-temperature PSD thus scales 
like S h (ui) oc ui 3 . This is in full agreement with the results 
obtained in Refs. 140114 II . where this asymptotic behavior 
is found to be uniform throughout the whole frequency re- 
gion. However, this is no longer so when r is finite: the 
second term in the rhs of Eq. ((33) introduces a linear cutoff 
in the limit ui — > 0, S h (u>) oc lj. In distinct contrast, in the 
high-frequency region, the first term in the rhs of Eq. (|33] 
becomes dominating. As a consequence, the PSD d33l ap- 
proaches a square-law asymptotic crossover dependence, 
o(oj) oc w 2 , in the high-frequency limit. 

4 Conclusions and open sundry topics 

By using the Green function formalism we have investi- 
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gated electronic heat current. Our focus centered on the is- 
sue of the heat current fluctuations in a molecular junction 
model composed of a single orbital molecular wire. For the 
noninteracting case we succeeded in deriving a closed form 
for the frequency dependence of heat current noise; i.e. the 
heat noise PSD, both in nonequilibrium XL ^ Tr and in 
thermal equilibrium Tl = Tr . The dependence of the heat 
current noise on the orbital energy £o is qualitatively differ- 
ent from that for the accompanying electric current noise, 
see Fig. [5] Moreover, the heat current fluctuation proper- 
ties depend strongly on the the overall tunneling coupling 
strength T L = T R = T. 

In the zero-temperature limit, the PSD of the heat cur- 
rent noise obeys two distinctive asymptotic behaviors, be- 
ing different in the low-frequency and the high-frequency 
regions. The particular square-law shape of the PSD in the 
high-frequency region is due to the Lorentzian shape of the 
transmission coefficient T{E) in Eq. (F29b . Yet, the gen- 
eral effect would remain for any choice of the coefficient 
in the form of a localized, bell-shaped function: the noise 
spectrum will deviate from a cubic power-law asymptotic 
behavior upon entering the high-frequency region. 

As emphasized in our introduction, with this work only 
the electron subsystem has been considered. Realistic heat 
transport in real molecular junctions would involve the 
complexity of interacting electrons and electron-phonon 
interactions [6|. This electronic heat transport may dom- 
inate in certain situations so that the measured heat noise 
can be attributed approximately to the electronic compo- 
nent only. The unified approach, which would include both 
the electron and the phonon subsystems, as well as the 
effects of their interactions, presents a future challenge 
although several contributions in this direction for the av- 
erage heat current (but not the heat current noise PSD) 
have already been undertaken before |i6. 25, 28]|48l . 

4.1 Open issues 

We conclude this study with further remarks that may shed 
light on pressing open problems and in addition may in- 
vigorate others to pursue future work along these lines. 
A first observation is that we obtained within the Green 
function analysis tractable expressions for quantum trans- 
port in the steady state without ever having to invoke the 
explicit knowledge of the inherent nonequilibrium density 
operator. Naturally, the quantum averages for the current 
and the auto-correlation of the quantum fluctuations carry 
less information as encompassed with the full steady state 
nonequilibrium density operator. The latter nonequilibrium 
density operator is typically very difficult to obtain and ex- 
plicit results are known for tailored situations only. In fact, 
explicit results are very intricate already for those cases 
with overall quadratic Hamiltonians only [4-9 1 . 

A much more subtle issue refers to the experimental de- 
tection of quantum correlation features. In clear contrast to 
the case with a quantum, single-time expectation of a quan- 
tum observable the issue of measurement of manifest quan- 



tum correlations is a delicate and difficult issue that is only 
rarely addressed with sufficient care in the literature. This 
is so because the mere calculation of a theoretical two-time 
quantum correlations does not say anything about its feasi- 
ble experimental measurement scenario. Either strong, i.e. 
of the von Neumann-type, or weak quantum measurements 
impact the dynamics as clearly manifested with the exam- 
ple of with the Zeno-effect |50 5 1 1. With more that one 
time present this objective relates to the problem of mea- 
surements of quantities that are not quantum observables 
[46 47 52 53 1. To appreciate the complexity somewhat in 
more detail let us first consider the case with classical ran- 
dom variables. Then the PSD can be obtained experimen- 
tally as the limit of a time average of the classical random 
process J'l(t), via considering the 

1 



S(uj) =lim T ^ 00 ^ 



exp(iut) dt 



(34) 



Note that classically the measurement of the random vari- 
able of the momentary heat flow J h (t) at time t presents 
no serious conceptional problem while the same is not at 
all straightforward as a function of continuous time t for a 
quantum dynamics. Moreover, even classically, the result 
in Eq. ( f34b holds true only when the finite value St{u) 
tends to the exact value S(uj) with its variance approaching 
zero as T — >• oo. The latter implies conditions of higher, 
fourth-order correlations to be satisfied [E>4| . With the non- 
commutation property of quantum observables at different 
times at work no such direct scenario is available for the ex- 
periment. Here the complexity of quantum measurements 
will enter in its full generality. In certain situations, how- 
ever, the task simplifies fortunately. 

Consider the case of quantum linear response theory. 
The measurement of a single observable (here the heat flux 
operator) due to an external perturbation is typically re- 
lated to the evaluation of a specific quantum correlation 
function. The case of the quantum-dissipation relation of 
Callen-Welton presents such a celebrated case Q55ll56ll57l . 
There, the dissipative part of the measurable, frequency- 
dependent susceptibility of a perturbed observable B is 
uniquely related to the PSD Sba(w) of quantum fluctua- 
tions of the observable B and the fluctuations of observable 
A to which the applied external force couples. In our case it 
remains therefore a formidable task to research for the cor- 
responding variable so that the single-time measurement 
of its linear response becomes related to the PSD in Eq. 
(|30| >. This at best seems possible for thermal equilibrium 
PSD in which an imposed energy perturbation couples to 
the thermal affinity AT/T. Seemingly this is not possible, 
however, for the equilibrium heat flow fluctuations at abso- 
lute T = with an inherent divergent affinity B2l . In pres- 
ence of quantum coherence destroying phenomena such as 
high temperatures or disorder the nature of quantum cor- 
relations become suppressed. Then, the classical scenario 
can be used again to verify the theoretical predictions ||3] 

E2. 
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